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ABSTRACT 


It  is  the  purpose  of  this  paper  to  develop  a  useful  mathematical 
model  of  ASW  aircraft  availability.  The  increasing  emphasis  of  systems 
studies  dictates  the  use  of  accurate  and  representative  models  of  the  ASW 
systems.  At  present,  many  studies  are  using  essentially  the  same  models 
developed  during  World  War  IL  This  paper  is  an  attempt  to  make  use  of 
advanced  theory  in  a  more  powerful  and  flexible  model  and  to  make  the  use 
of  the  model  practical  and  verifiable. 

The  writer  adapted  the  time  homogeneous  bivariate  model  as  de¬ 
veloped  by  F.  C.  Collins.  This  is  a  discrete  time  Markov  process  with 
a  stochastic  matrix  of  transition  probabilities  wherein  the  maintenance 
process  is  modeled  ar  a  pulsed  input  multiple  server  queue. 

The  model  was  programmed  in  FORTRAN  63  on  the  CDC  1 604  and 
then  modified  to  illow  for  variability,  in  the  input  parameters.  Other 
modifications  include  an  increase  in  the  size  of  the  model  to  accommodate 
a  1 6-aircraft  squadron,  the  largest  ASW  squadron  at  present,  and  an 
explicit  form  solution  to  the  maintenance  queueing  equations. 
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1.  INTRODUCTION 


The  threat  to  freedom  of  the  seas  posed  by  the  vast  Soviet  sub¬ 
marine  fleet  is  perhaps  the  most  thorny  problem  facing  the  U.  S.  Navy 
today.  Two  world  wars  have  produced  Pyrrhic  victories  over  limited 
submarine  fleets.  During  the  Second  World  War  operations  analysis 
was  born  into  the  Navy  to  aid  in  the  defeat  of  the  German  submarine. 

The  classic  antisubmarine  warfare  (ASW)  analyses  and  models  developed 
by  Morse  [2]  and  Koopmans  [3]  are  still  being  used  today,  over  two 
decades  later,  in  most  of  the  ASW  study  efforts  for  the  Navy. 

These  early  ASW  analyses  assumed  a  given  level  of  search  effort 
available  and  directly  evaluated  the  probability  that  an  ASW  subsystem 
could  detect  and/or  kill  a  submarine.  This  assumption  is  not  only 
logical  to  make  the  problem  tractable,  but  also  practical  since  no 
immediate  changes  in  ASW  force  levels  could  be  expected.  Moreover, 
the  studies  were  conducted  during  the  war,  not  before  it  started.  It  is 
the  purpose  of  this  paper  to  present  a  probabilistic  model  to  describe 
the  available  effort.  Such  a  model  can  be  used  to  sharpen  the  estimates 
of  the  effectiveness  of  an  ASW  subsystem  and  to  study  the  charac¬ 
teristics  of  the  associated  support  system. 

Naturally,  the  current  study  plays  an  important  but  limited  role  in 
the  overall  problem  of  designing  an  entire  ASW  system.  The  diffi¬ 
culties  involved  in  such  a  specification  are  legion.  First  and  foremost 
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is  the  quantification  of  the  ASW  mission  in  denying  the  enemy  the 
effective  use  of  his  submarines.  Currently,  the  probability  of  detecting 
and/or  killing  submarines  is  used  as  the  measure  of  effectiveness  of 
the  mission,  and  it  appears  that  a  more  encompassing  one  has  not  been 
developed.  Second,  the  specification  of  an  ASW  force  level  to  counter 
a  given  threat  has  many  inherent  subjective  elements.  These  are  due  to 
the  existing  historical  bias  in  predicting  the  conduct  of  a  future  ASW  war 
with  an  enemy,  particularly  one  who  has  never  before  used  a  large 
submarine  force  in  its  military  operations.  The  reader  can  imagine 
why  merely  defining  terms  such  as  •'threat11  and  "effective  counter" 
becomes  quite  difficult. 

Thus,  there  is  a  need  to  investigate  the  levels  of  search  effort 
specified.  This  may  require  acceptable  models  to  measure  the  availa¬ 
bility  of  effort,  its  effectiveness,  and  determine  the  logistic  support 
required  for  any  level  of  available  effort.  Specifically,  the  ASW  sub¬ 
system  to  be  modeled  is  the  carrier -based  aircraft,  although  the  model 
is  adaptable  to  other  systems. 

The  method  of  investigating  the  demand  for  ASW  carrier  a/c  will 
assume  that  the  desired  number  of  a/c  on  station  is  known  as  an  input 
parameter.  The  support  required  to  achieve  this  measure  of  available 
effort  depends  upon  maintenance  space,  manpower,  and  supply. 
Generally,  we  shall  consider  how  an  ASW  carrier  supports  this  number 
of  a/c  on  station  with  the  present  or  proposed  number  of  a/c  embarked 
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on  the  carrier.  The  parametric  input  can  be  subjected  to  sensitivity- 
analyses. 

The  operational  commander  of  the  ASW  force  launches  the  desired 
number  of  a/c  on  station  to  screen,  search,  or  actively  prosecute  a 
submarine  contact.  Each  a/c  is  relieved  on  station.  Each  such  relief 
requires  the  launching  of  another  a/c  prior  to  the  recovery  of  the  initial 
a/c.  The  returning  a/c  must  receive  varying  degrees  of  maintenance 
and  requires  refueling  and  rearming.  This  cycle  continues  until  the 
mission  is  completed.  Loss  of  a/c  due  to  accidents,  insufficient  supply, 
and  lack  of  repair  capability  cause  deviations  in  this  procedure.  Naval 
operations  involve  the  interaction  of  many  quantities  which  are  random 
in  nature.  Not  all  can  be  considered  in  a  tractable  mathematical  model. 
Some  quantities  which  are  important  are  omitted.  One  example  is  the 
length  of  each  cycle  time,  which  is  assumed  to  be  a  constant  value. 
Including  variables  of  this  nature  incurs  unnecessary  mathematical 
complication.  It  is  hoped  that  adequacy  of  the  model  can  be  measured 
by  using  fleet  data  available  from  the  Fleet  ASW  Data  Analysis  Program 
(FADAP). 

Collins  [5]  describes  a  bivariate  Markov  model  for  airborne  early 
warning  (AEW)  and  combat  air  patrol  (CAP)  jet  a/c  operating  in  an 
attack  carrier  force.  This  model  is  used  to  evaluate  the  probability  of 
maintaining  a  fixed  requirement  of  a/c  on  station  as  a  measure  of  ef¬ 
fectiveness  of  the  system.  It  has  subsequently  been  used  in  a  larger 
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attack  force  study  for  the  Navy.  The  model  computes  the  probabilities 
of  the  number  of  a/c  on  station  and  in  or  awaiting  maintenance  at  any 
given  launch  period.  The  comparable  ASW  problem  differs  in  the  following 
aspects: 

1.  Type,  range,  and  speed  of  a/c; 

2.  The  variable  number  of  a/c  required  for  mission; 

3.  Attrition  due  to  accidents  and  supply  failures; 

4.  The  greater  number  of  ASW  a/c. 

It  was  decided  to  use  the  Collins'  model  with  appropriate  modification. 

For  immediate  reference,  the  mathematical  content  of  the  model  will 
be  repeated  herein. 

In  order  to  incorporate  these  modifications,  it  was  necessary  to 
spend  some  time  reprogramming  on  the  COC  1604  digital  computer  in 
FORTRAN  63,  the  CDC  version  of  the  IBM  FORTRAN  IV.  The  original 
program  was  not  readily  available  and  was  written  in  an  early  assembler 
language.  Moreover,  the  numerical  analysis  was  not  sufficiently  sharp 
to  handle  the  larger  input  values.  Also,  double  precision  (two  computer 
words  instead  of  one)  arithmetic  was  required  in  one  subroutine  for  an 
accurate  explicit  solution  to  the  maintenance  queueing  equations  (see 
Appendix  I).  This  effected  a  50%  decrease  in  the  computer  time 
required  for  developing  a  matrix  of  transition  probabilities. 

Following  this  introduction,  section  2  contains  a  brief  description  of 
the  operational  problems  involved  and  the  assumptions  made.  A  brief 
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description  of  Markov  chains  and  the  mathematical  model  are  presented 
in  section  3.  The  details  for  computing  the  matrix  of  transition  proba¬ 
bilities  are  given  in  section  4.  General  employment  of  the  model 
follows.  The  appendices  include  the  solution  mentioned  on  the  pre¬ 
ceding  page,  a  logical  flow  diagram  of  the  program,  a  copy  of  the 
program,  and  some  sample  results. 
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ASSUMPTIONS 


The  real-world  employment  of  carrier  a/c  is  cyclic  in  nature,  and 
the  present  state  of  any  given  a/c  (i.  e.  ,  flying,  in  or  awaiting  mainte¬ 
nance)  depends  largely  on  what  the  previous  state  was.  This  fact 
suggests  that  a  Markovian  assumption  can  logically  be  made  for  the 
a/c  transition  probabilities.  In  the  search  phase,  a/c  may  or  may  not 
relieve  on  station;  but,  in  any  part  of  the  contact  investigation  phase, 
relief  on  station  will  be  made.  To  insure  full  screening  and  mission 
coverage,  a/c  will  relieve  on  station. 

The  question  of  resupply  during  an  operation  depends  primarily  on 
the  availability  of  carrier  on-board  delivery  (COD).  This  depends  on 
the  geographical  location  and  the  mission  (convoy  protection,  strike- 
force  protection,  huncer -killer  operation,  etc.).  In  practice,  resupply 
is  not  anticipated  within  a  week's  period,  and  around-the-clock  oper¬ 
ations  have  continued  for  two  weeks  without  resupply. 

Standard  maintenance  procedures  aboard  carriers  preclude  major 
maintenance  on  the  flight  deck.  It  will  be  assumed  that  sufficient  notice 
is  given  so  that  all  major  120-hour  checks  will  be  completed  prior  to 
the  operation.  This  assumption  can  be  modified  with  an  appropriate 
adjustment  in  the  mean  repair  rate.  The  concept  of  maintenance  crews 
assigned  to  hangar  deck  areas  ("spots"),  as  developed  by  Collins  [  3 j , 
will  be  used.  Each  crew  will  be  capable  of  all  types  of  maintenance 


and  will  operate  independently  at  the  identical  mean  repair  rate  X. 

The  number  of  spots  is  determined  by  the  average  number  of  such 
crews  available  to  work  continuously  around  the  clock  on  a  watch  basis. 

The  state  of  each  a/c  is  assumed  to  be  statistically  independent  of 
that  of  others,  and  the  launching  and  landing  transition  probabilities 
will  be  developed  on  the  basis  of  independent  Bernoulli  trials.  The 
parameters  can  be  determined  using  the  maximum  likelihood  estimators. 
The  range  of  the  number  of  a/c  desired  on  station  at  any  given  cycle  will 
be  set  by  the  user.  The  number  to  be  launched  at  any  time  is  assumed 
equally  likely  within  this  range.  This  input  parameter  is  a  function  of 
the  estimated  submarine  density  (i.  e.  ,  expected  contact  rate).  The 
lower  limit  will  be  set  at  the  number  of  a/c  desired  on  station  in  the 
search  (screening)  phase,  and  the  upper  limit  is  set  at  the  maximum 
practicable  number  of  a/c  to  be  launched  during  a  multiple -contact 
phase. 

Briefly,  the  assumptions  are: 

1.  a/c  will  be  relieved  on  station. 

2.  Any  desired  length  of  operation  can  be  set  as  an  input. 

3.  Major  120 -hour  checks  will  be  completed  prior  to  the  operation. 

4.  No  resupply  to  th?  carrier  is  available. 

5.  The  launch-to-launch  cycle  for  all  ASW  a/c  is  four  hours, 

6.  Minor  maintenance,  refueling,  and  rearming  only  can  be 
performed  on  the  flight  deck. 
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7.  Each  maintenance  spot  is  characterized  with  an  independent 
exponential  repair  time  with  mean  repair  rate  of  X  for 
around-the-clock  operations. 

8.  The  number  of  a/c  lost  due  to  attrition  is  a  Poisson  random 
variable  for  each  cycle  period  with  parameter  X  (a/c 

A 

accident/flying  hours  for  a/c  type). 

9.  Any  a/c  lost  by  accident  will  not  be  returned  to  service  due  to 
either  (a)  physical  loss  at  sea,  or  (b)  insufficient  mainte¬ 
nance  capability  aboard  ship  and  lack  of  major  parts. 

10.  The  number  of  a/c  launched  for  each  cycle  is  uniformly 

distributed  between  the  upper  and  lower  limits  determined  by 
the  user. 
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3. 


MODEL  DESCRIPTION 


3.  1  The  Theory 

A  stochastic  or  random  process  is  a  collection  of  random  variables 
indexed  on  some  set  T,  (X(t)  ,  t«  T).  In  this  case,  time  is  the  indexing 
set,  and  the  Markovian  assumption  states  that  the  future  state  of  the 
process  depends  only  on  the  state  at  the  present  time  and  not  on  its 
past  history.  Due  to  the  cyclic  nature  of  our  problem,  it  is  possible 
to  increment  time  (T  =  (0,  1,  . .  .  ) )  using  the  cycle  time  from  launch 
to  launch  as  the  steps  of  unit  time  in  a  discrete  Markov  chain.  It  is 
assumed  that  the  reader  is  familiar  with  the  notion  of  a  random  variable 
as  a  function  defined  on  a  sample  description  space  (S)  on  which  the 
family  of  events  or  outcomes  (E)  of  a  probability  function  can  be 
defined  [4] . 

A  discrete  time  Markov  chain  is  described  by  a  sequence  of  dis¬ 
crete  valued  random  variables  and  is  determined  when  the  one -step 
transition  probabilities  of  the  state  variables  are  specified,  i.  e.  ,  a 
conditional  transition  probability  of  a  transition  at  time  n  for  each  pair 

of  i,  j  -  0,  1 . m  (m  being  the  number  of  states  in  the  process) 

must  be  given. 

p..(n,  n  +  1)  =  P  [X  (n  +  1)  =  j  |  X(n)  =  i] 

If  the  transition  probability  functions  depend  only  on  the  time  difference, 
we  have  time  homogen  uty 
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p^(n  +1,1)  =  p„(0,  1)  =  • 

The  initial  state  of  the  system  must  be  given  either  as  a  specific  state 
or  randomly  as  a  probability  distribution  function  over  the  possible 
states. 

The  p..  (transition  probabilities)  are  arranged  in  matrix  form  and 
satisfy: 


1. 

2. 


p.  .  2  0 

1J 


m 

£ 

j=0 
for  all 


PU 


for  i,  j  =  0,  1 . m  ; 

=  1  ,  i,  e.  ,  the  rows  of  the  transition  matrix  sum  to  1 
i  for  the  states  within  the  description  space  [4]. 
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3.  2  The  Model 


In  order  to  establish  the  finite  set  of  states  (E)  for  the  model,  we 
shall  consider  two  random  variables  defined  as  follows: 

Xj  (t)  =  The  number  of  a/c  flying  at  time  t  not  having  flown  in 
the  previous  launch -to -launch  interval. 

X£(t)  =  The  number  of  a/c  in  or  awaiting  maintenance  at  timet. 

Now,  we  will  consider  the  vector  X(t)  =  [X^  (t) ,  X^  (t)  ]  as  a  pair  of 
random  variables  and  thereby  have  a  bivariate  stochastic  process  with 

*  3 

the  possible  states  ranging  from  (0,  0)  to  (A,  N). 

0  s  Xj  (t)  s  A  =  No.  of  a/c  desired  on  station,  and 

0  £  X^ft)  £  N  =  No.  of  a/c  of  given  type  aboard  carrier. 

We  will  define  an  operating  cycle  as  an  interval  unit  of  time.  Process 
observations  of  X(t)  will  be  made  at  successive  unit  interval  launch 
times.  To  develop  the  p„  elements,  consider  a  given  time  t  for  launching 
until  A  aircraft  are  flying  or  until  the  supply  of  r^ady  a/c  is  depleted. 
Those  a/ c  failing  the  launch  enter  the  maintenance  state  at  this  idealized 
point  in  time  t  (the  total  launching  time  required  is  much  less  than  the 
total  cycle  time)i  At  some  time  T,  less  than  the  launch-to-launch  unit 
time  interval,  the  a/c  which  were  relieved  on  station  return  and  land  at 
the  idealized  point  in  time  t  +  T.  Some  of  these  a/c  will  require  mainte¬ 
nance  and  enter  the  maintenance  queue.  Those  requiring  only  refueling 
and  preflight  inspection  will  enter  a  ready  status  to  be  tested  for  the  next 
launch. 
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During  the  unit  time  interval,  maintenance  will  be  performed  on 

those  a/ c  in  the  not-ready  status,  and  a  certain  number  of  aircraft  will 

be  repaired  according  to  assumption  7. 

In  summary,  we  start  the  system  in  some  initial  state  (such  as 

(0,  0)  with  no  a/c  flying  or  in  maintenance)  or  start  with  a  probability 

distribution  Q(t  )  over  the  states,  E,  at  time  t  .  We  launch,  recover, 
o  o 

and  repair  a/c  in  the  unit  interval  and  repeat  the  process  over  each 

succeeding  unit  time  interval  until  the  end  of  the  operating  period. 

Knowing  the  transition  probabilities  within  the  unit  time  interval,  we  can 

develop  the  elements  of  the  transition  matrix,  P,  or  {p.  ..  .  . 

i  l)  i  \  P»  J; 

These  are  the  probabilities  of  going  from  the  state  of  a  a/c  flying  and 
i  a/c  in  maintenance  to  £  a/c  flying  and  j  a/c  in  maintenance  over  the 
unit  time  interval. 

It  was  assumed  in  section  2  that  A,  the  number  of  a/ c  to  be  launched, 
and  N,  the  total  number  of  a/c  on  board,  are  random  variables,  whereas 
they  have  been  treated  as  constants  so  far  in  the  development.  To  be 
analytically  correct  in  including  this  feature,  one  should  develop  the 
appropriate  quadrivariate  process.  Such  a  development  leads  to  too 
large  a  state  space  and  the  author  chose  to  include  these  effects  by  using 
a  Monte  Carlo  simulation  technique.  That  is,  at  the  beginning  of  each 
cycle,  a  random  mechanism  is  used  to  determine  the  values  on  A  and  N. 

The  probability  of  iesing  an  a/c  or  changing  the  desired  number  to 
be  launched  is  determined  from  the  specified  distributions  at  the  beginning 
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of  each  unit  interval,  and  the  resulting  P  matrix  containing  the 

p  t  is  then  recomputed.  The  probability  distribution  Q(t) 

(a,  1).  (P  ,  j) 

over  the  states  at  any  time  t  may  be  determined  by  the  appropriate 
number  of  successive  iterations  of  the  Q  vector  times  the  P  matrix,  i.  e.  , 
Q(t)  =  P[XL(t)  =  p,  X2(t)  =  j]  =  Q(t  -  l)xP  . 

The  probability  of  maintaining  a  a/c  on  station  over  any  given  period 
of  operation  may  be  obtained  at  any  unit  time  t  (i.  e.  ,  the  beginning  of 
the  next  cycle)  by  summing  out  the  appropriate  maintenance  state  proba¬ 
bilities.  Thus,  P(of  a/c  are  flying  at  time  t  )  = 

Pr  (X  (t)  *  or  )  =  £  Pr  (X  (t)  =  a ,  X2(t)  =  i  )  . 
i=o 

A  mathematical  comment  appears  to  be  in  order.  In  the  case  of 
fixed  A  and  N,  the  states  of  the  Markov  chain  are  positive  recurrent;  and 
steady- state  probabilities  can  bo  found  for  the  entire  state  space.  In  the 
case  of  decreasing  N  due  to  a/c  attrition,  this  13  not  true;  and  (0,  0) 
becomes  an  absorbing  state  as  time  (t)  goes  to  infinity.  This  latter 
consideration  is  not  a  realistic  one  for  the  operational  period  envisioned. 
Therefore,  it  is  mathematically  more  feasible  to  use  the  former  chain  in 
conjunction  with  the  Monte  Carlo  technique. 
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4.  DEVELOPMENT  OF  THE  TRANSITION  MATRIX 


Perhaps  the  simplest  way  to  view  this  developmer,  is  to  note  the 
various  transition  probabilities  incorporated  in  one -unit  time  cycle 
defined  as  follows: 

(1)  y  =  the  launching  transition  probabilities  at  time  t.  This  is  the 
fgh 

probability  of  taking  f  ready  a/c,  launching  g  successfully,  and  sending 
h  into  maintenance.  Each  a/c  to  be  launched  is  considered  a 
Bernoulli  trial  with  probability  of  failure  of  p  ,  which  is  estimable 
and  subject  to  sensitivity  analysis.  The  values  of  are: 

a.  0  if  g  >  A,  since  only  A  a/c  are  desired; 

b.  0  if  g  +  h  >  f ;  it  is  impossible  to  launch  and  send  into  mainte¬ 
nance  more  a/c  than  are  available; 

c.  0  if  g  <  A,  g  +  h  <  f ;  launching  continues  until  A  a/c  are  flying 

or  until  all  f  are  used  up; 
f  o  f  -  c 

d.  (  )  ( 1  -  p  )  (p  )  6  if  g  <  A,  g  +  h  =  f ,  standard  binomial 

®  Y  Y 

when  all  a/c  in  the  ready  state  are  used  up  but  the  A  a/c  are 

not  launched; 

2  +  h  -  1  b  h 

e.  (  k  )  (1  -  p)  6  (p)  if  g  =  A,  g  +  h>  f ,  standard  negative 
binomial  for  g  successes  in  g  +  h  -  1  trials. 


{2)  n<*  (m)  =  the  landing  transition  probabilities  which  occur  at  time 

t+  T  .  We  must  consider  the  probability  that  if  there  ar i  a/c  flying 
at  time  t  then  m  a/c  will  enter  maintenance  at  recovery  time  t  +  T . 
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II  (m)  will  equal  a  standard  binomial  where  p  =  the  probability  of 
a 

equipment  failure  in  flight: 

_  .a  .  . o' -m,  m  _  , 

II  (m)  =  (  )  (1  -  p)  (p)  ,  m  =  0,  1,  .  . .  ,  a  . 

a  m 

p. .  (t  )  =  the  maintenance  transition  probabilities,  i.  *. ,  the  proba- 

bility  of  repairing  (i  -  j)  a/c  in  time  t  .  Two  maintenance  periods 

occur:  the  first  starting  at  time  t  and  ending  at  time  t  +  T  ,  the 

second  starting  at.  time  t  +  T  and  ending  at  the  end  of  the  cycle, 

(t  +  1).  Under  assumption  7,  the  pulsed  input,  multiple  exponential 

server  queue  is  developed  with  D  maintenance  "spots"  or  servers 

each  with  identical,  independent  service  rates,  X  .  For  each  server, 

then,  the  probability  of  remaining  occupied  (given  the  server  is 

\  T 

busy)  in  time  t  =  e  .  The  probability  of  becoming  free  (i.  e.  , 

“  X  T 

repairing  an  a/c)  =  1  -  e  .  The  resulting  queueing  equations  are: 

A.  dP.  (t)  /  dt  =  -  n  X  P.  (t)  +  (n+l)XP.  ^  ,  (t)  forO£n<D; 

i,  n  l,  n  l,  n+  I 

B.  dP.  (t)  /  dt  =  -  DX P.  (t)  +  DXP.  Ll(t)  for  n  2  D  . 

l,  n  i,  n  l ,  n+  i 

Three  ranges  of  i  (initial  queue  state) ,  j  (final  queue  state) ,  and 
D  become  significant: 

a.  When  j  £  i  £  D,  then  not  all  spots  are  busy  since  there  are 
fewer  a/c  in  maintenance  than  spots.  Each  spot  works  inde¬ 
pendently;  therefore,  the  solution  to  A  is  the  binomial: 


b.  When  Ds  j  S  i,  then  all  spots  are  occupied  throughout  the  total 
service  time,  and  the  closed  form  solution  to  B  is  the  Poisson: 


c,  When  j  <  D<  i,  then  all  spots  are  busy  at  the  beginning  of  the 
service  period,  and  some  spots  become  idle  during  the  service 
period.  The  explicit  form  solution  of  equation  A  is  found  using 
moment  generating  function  transformation: 


Py(t) 


*  <"><?>  {(ofr) 


(i  -  D)  -  \t  n 


*  (^r)k(5^) 


i  -  D 


(The  derivation  of  this  solution  is  discussed  in  Appendix  I. ) 


The  figure  on  the  following  page  will  show  the  relationships  of  these 
transition  probabilities  within  the  unit  time  interval. 
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TRANSITION  PROBABILITIES  WITHIN  THE  UNIT  CYCLE 


FIGURE  1 


In  order  to  develop  each  transition  probability  over  the  total  unit 
time  interval,  we  must  consider  all  events  taking  place  within  the 
interval.  Thus,  to  obtain  the  probability  of  going  from  a  a/c  flying  and 
i  a/c  in  maintenance  to  3  a/c  flying  and  j  a/c  in  maintenance,  we  start 
at  the  state  ( a ,  i)  at  time  t.  At  this  time,  a/c  are  launched  and  some 
1  a/c  failing  the  launch  enter  maintenance.  These  i  +  1  in  maintenance 
are  then  serviced  until  time  t  +  T  when  some  k  a/c  are  still  in  the 
maintenance  state.  At  time  t  +  T,  of  the  a  a/c  previously  flying,  some 
m  enter  maintenance  and  (a  -  m)  enter  the  ready  pool.  Maintenance  is 
continued  on  the  (k  +  m)  a/c  for  the  remainder  of  the  cycle  (1  -  T)  ,  until 
the  end  of  the  unit  time  interval  when  j  a/c  remain  in  the  maintenance 
state.  In  functional  form: 


(a,  i) ,  (3  ,  j) 


N-o-i  i+1 

a 

r  l 

L 

l=o  k=o 

m=o 

•P 

i+  1  ,  k 

N-a-i,  3,  1 


P  x  i  v.  (T)  •  n  (m)  .  p,  x  .  (1  - 
i+  1  ,  k  a  k  +  m  ,  j 


T) 
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5.  SUMMARY 


Representative  values  for  the  mean  repair  rate  and  the  landing  and 
launching  failure  rates  produced  results  in  agreement  with  the  sensitivity 
analysis  by  Collins  on  these  parameters  in  [5],  For  failure  proba¬ 
bilities  less  than  .  5,  and  mean  repair  rate  less  than  12  hours,  the  effect 
of  reducing  the  available  maintenance  time  to  80%  of  the  cycle  time  was 
negligible.  Optimal  loading  and  cycling  policies  can  be  determined  for 
known  values  of  these  rates. 

The  model  affords  the  following  checks:  (1)  the  rows  of  each  P 

matrix  are  summed  as  they  are  computed  by  the  program;  and  (2)  the 

probability  distribution  vector  (Q  J)  is  summed  over  the  states.  Each 

-  8 

summation  was  within  10  of  one  in  the  computer  model. 

* 

The  user  may  substitute  any  available  distribution  over  the  interval 

of  a/c  desired  on  station.  In  order  to  keep  A  fixed,  enter  the  desired 

value  as  both  upper  and  lower  limit  ( A  =  ALOLIM  =  LUPLIM) .  For 

-  8 

fixed  N,  use  a  very  small  value  for  ALAM  (such  as  10  ).  Subroutine 

KRAN  is  a  uniform  generator,  using  the  half  open  interval  (lower  limit  +  1, 
upper  limit  +  2)  and  a  starting  number  as  inputs.  KRAN  outputs  an  integer 
in  this  interval.  Subroutine  DRAW  was  used  to  provide  some  intuitive 
grasp  of  the  results.  DRAW  was  used  in  binary  card  form  and  is  not 
essential  to  the  main  program.  (The  indicated  associated  statements 
must  be  removed,  however. ) 
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The  results  of  reasonable  arbitrary  parameter  values,  based  on  the 
author's  experience,  have  shown  that  most  of  the  probabilities  concen¬ 
trate  over  a  few  states.  Moreover,  computation  time  increases  rapidly 
as  a  function  of  N  (no.  of  a/c),  see  Figure  2.  This  would  indicate  that  a 
simple  approximation  to  the  model  could  be  developed.  One  method 
presently  being  investigated  to  reduce  comput?tion  time  is  to  shrink  the 
probability  state  space  to  include  only  those  significant  states  and,  thus, 
reduce  the  size  of  the  transition  matrix.  Alternatively,  the  eigenvector, 
eigenvalue  representation  of  the  P  matrix,  might  be  used. 

Originally,  it  was  hoped  to  utilize  the  data  from  the  Fleet  ASW  Data 
Analysis  Program  (FADAP)  to  attempt  a  verification  of  the  model  with 
its  real-world  counterpart.  The  only  method  available  at  present  for 
obtaining  the  necessary  data  is  by  direct  observation  or  a  program  of 
data  collection,  as  suggested  by  Collins  [  5]  . 

Many  fruitful  areas  of  investigation  exist: 

(1)  Attrition  has  been  simply  modeled  by  the  Poisson  method.  The  two 
components  of  attrition,  accidents  and  supply  shortage,  can  be  more 
accurately  modeled  and  used  to  develop  logistic  schedules  for 
maintenance  and  supply.  One  simple  technique  is  to  assume  each 
component  is  independent  and  Poisson,  and  estimate  a  supply  failure 
rate  for  AOCP  attrition  from  past  data.  With  these  assumptions, 
the  total  attrition  is  Poisson,  with  the  parameter  equal  to  the  sum 
of  the  accident  and  supply  failure  rates. 


30 


(2)  The  model  can  be  modified  to  make  the  number  of  maintenance  spots 
available  fcr  any  cycle  a  variable  function  of  time,  D(t) , 

(3)  An  investigation  of  the  Markovian  assumption  validity  as  the  cycle 
times  become  smaller  and  smaller. 

(4)  Development  of  a  continuous  time  model. 

(5)  Modification  of  the  model  to  simulate  resupply  by  COD. 

(6)  A  study  of  the  distribution  of  submarine  contacts  to  determine  the 
validity  of  the  uniform  a/ c  demand  assumption. 
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m 


in 


Time  (min) 


PROGRAM  ASSEMBLY  AND  COMPUTATION  TIME 
FOR  ONE  TRANSITION  MATRIX  (P)  AS  A  FUNCTION 
OF  THE  TOTAL  NUMBER  OF  AIRCRAFT  (N) 

FIGURE  2 
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APPENDIX  I 


EXPLICIT  SOLUTIONS  OF  THE 
MAINTENANCE  QUEUEING  EQUATIONS 

The  queueing  equations  for  the  pulsed  input  queue  are  essentially 
the  pure  death  process  given  in  [  1  ]  and  [4]  as  problems  and  developed 
by  Collins  in  [5].  The  equations  are: 


A. 


dP.  (t) 

i.  R 

dt 


-  n  A. P.  (t)  +  (n  +  1)  AP.  ,  ,  (t)  for  0  <  n  <  D 
x,  n  l,  n  +  I 


B. 


dP  (t) 
i»  n 

dt 


-  n  A.  P.  (t)  +  DAP.  ,  (t) 
l,  n  x,  n  +  1 


for  n  ^  D 


where  P..(0)  =  A.,  and  P..(t)  =  0  for  i<  j,  since  no  input  (arrivals) 
xj  xj  xj 

occur  during  the  service  time. 

Equation  B  is  solved  directly  in  closed  form: 


Pi  n(t) 

x,  n 


(DAt)  (n  "  X)  e  "  X  Dt 
(n  -  i)  ! 


Now  transforming  the  first  equation  (A)  using  the  moment  gener¬ 
ating  function  (MGF),  f 

D-l 

G(s,  t)  =  £  s  P  (t)  , 

«  n 

n=0 

as  outlined  in  [4]  (Chapter  7),  and  its  partial  derivatives: 
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(1) 


dG 


DL  sn  p'  (t) 
dt  n=0 


(2) 


dG  DZ}  n  -  1  ... 

=  Z;  ns  P  (t) 

ds  n 

n=0 


Where  (t)  denotes  the  conditional  probability  ^(t)  ,  by  substituting 
(A)  into  (1),  properly  identifying  the  first  summation  with  (2),  and 
changing  the  second  summation  index  to  r  =  n  +  1  ,  we  get: 


dG  .  dG  A  J  r  -  1  _  ... 

— T—  =  -  X  s  — —  +  K  Z  T  s  P  (t)  ,  or 

dt  ds  _  r 

r=0 


(3) 


dG 

dt 


-  X  (s  -  1)  ^  +  XDs  D  '  1  PD  (t)  , 


since 


D 

Z  r  s 
r=0 


r'lpr(t)  +  DS  °  ’  1  PD(t> 


Next,  replace  the  partial  differential  equation  (3)  with  a  system  of 

i" 

ordinary  differential  equations  using  the  Lagrangian  auxiliary  equations: 


ds 


dt 


1  \  (s  -  1) 


dz 


XDsD  ‘  1  PD(t) 


The  solution  to  the  first  equation  (using  the  first  two  differentials) 


is: 


Xt  =  In  (s  -  1)  +  C' 
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and  hence 


Xt 

8  =  C!e  +  1 
or 

Gj  =  e  '  U  (s  -  1)  . 

The  second  equation  is:  (using  first  and  third  differentials) 
dz  =  -\D  (C^e  Xt  +  1)D  '  1  PD(t)  dt  . 
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To  get  our  particular  solution,  use  the  boundary  conditions  for  G(s,  t)  : 


(1)  for  s  =  1  , 


G(  1 ,  t)  =  E  Pn(t) 
-  A 
n=0 


where 


=  Pr  [no.  in  maintenance  at  t  is  <  D  I  i  at  t  =  0] 


i-D  e"XDt(XDt\n 

G(l,  t)  =  1  -  E  - ~f— =  1  -  tj(t) 

n=0 


u  ( 1 ,  t,  z)  =  =  0 


f.  A  %  ,  -\Dt  ®  (XDt)k 

v  (i,  t,  z)  =  C  =  z  +  e  E  - ^ r - 

k=0 


C2  =  z  +  (*) 


(2)  for  t  =  0  , 


where 


G (s,  0)  =  £  s  P  (0)  =  0  ,  since  iz  n>  D 

~  n 

n=0 


u  (s,  0,  z)  =  (s  -  1) 


D-l  n  m 

v(s,  0,  z)  =  C2  =  Z  +  £  (  j  )  ( 3  -  1  )J  Q  -  ■ 


Thus, 


D-l  _  m 

v  /  n  j  / _ P  *\ 


G  (s,  0)  =  z  +  S  (?)  CJ  (^Tj) 


Substituting  the  general  value  for  above : 


D-l  _  .  ^  .  m 

G  (s,  t)  =  0  (u,  v)  =  L  (  7)  (s  -  1)J  e  ’ 

j=0 


(dTj) 


-°L  <?><.-  1)J  E  e-XDt 

j=0  J  k=0  K 


CdTj) 


m  -  k 


Rearranging  terms, 


D-l  _  D-l 


g(s,  t)  =  l  sn  l  0)  ( 7  )  <-l)J  f  ( -E-A 

n=0  j=n  n  3  L^D-jJ 


D  e-Xtj 


-  e 


m  -  k.' 


k  Dt  ™  (XDt)  z'  D  \m  “  -1 

J0  J 


where  P  (t)  =  the  coefficient  of  s 
n 
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COMPUTE  p{  >y)(t)  »  PR(I,XBETA,IY) 

(  I  -  IX+( IALPH-1 ) xTOTAL  A/C) 

f 


PRINT  THE  p±1*S  OF  THE  (IALP::,IX)TH  ROV.  AND  ROw  SUM 


APPENDIX  III 


THE  COMPUTER  PROGRAM 
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nnnnnnnnnnnnan nnnnnn 


-COOP » » LANMAN ,0 /49/S/ IS/ 2S/E/ 45=54 » 1  5  ,30000*5.  0001 

-BINARY, 56.  0001 

(REL0C0M.  0002 

-FTN,E.  0003. 

PROGRAM  MARKOV  0004. 

THIS  PROGRAM  IS  A  NONSTATIONARY  BIVARIATE  MARKOV  CHAIN  MODEL  OF  ASW  A/C  0005; 
OPERATIONS.  THE  RANDOM  VARIABLES  ARE  THE  NUMBER  OF  A/C  FLYING  AT  THE  0006! 
BEGINNING  OF  ANY  GIVEN  LAUNCH  CYCLE.  THE  MAXIMUM  NO.  OF  A/C  ALLOWED  IN  0007| 
THE  MODEL  IS  16<NA>.  THE  RANGE  OF  A/C  TO  BE  LAUNCHED  AT  ANY  GIVEN  0008. 

INTERVAL  IS  0  TO  6  A/C.  THE  FOLLOWING  INPUTS  ARE  REQUIRED.  0009 

I D=  THE  NO.  OF  INDEPENDENT  MAINTENANCE  SPOTS  0010 

NA=  TOTAL  NO.  OF  A/C  TYPE  ON  BOARD  0011 


TIME=TIME  FROM  LAUNCH  TO  RFCOVERY/ LAUNCH  TO  LAUNCH  CYCLE  TIMEf  HRS ) 00 12 ! 
FLAM=MEAN  REPAIR  TIME  PER  SPOT/LAUNCH  TO  LAUNCH  CYCLE  TIME  (HRS)0013‘ 
PGAM=  PROBABILITY  OF  A/C  FAILING  LAUNCH (M.L.EST.  FROM  PAST  DATA)  0014; 
p*  PROBABILITY  OF  A/C  FAILURE  DURING  FLIGHT  REQUIRING  MAINTENANCE  0015 


AT  LANDING  (M.L.  ESTIMATOR  FROM  PAST  DATA)  0016 

QIx  THE  PROBABILITY  DISTRIBUTION  VECTOR  OVER  ALL  POSSIBLE  STATES  0017. 
(7  X  ) 7  =  119)  SUCH  THAT  THE  SUM  OF  ALL  QI (I >  *  1.  THIS  0018 

IS  ESTIMATED  BY  THE  USER  AND  INPUTTED  BY  USING  A  DATA  STATEMENT0019 ■ 
ICYCLF  =  NO.  CYCLES  DESIRED  FOR  OPERATION  0020. 

v, CYCLE  x  LAUNCH  TO  LAUNCH  TIME(HRS)  <  TOT.  TIMExICYCLE  X  JCYCLE)  0021  ; 
ALAM  =  ACCIDENT  RATE  FOR  TYPE  A/C  (ACCIDENT/HOURS)  0022 

ALOLIM  =  DESIRED  LOWER  LIMIT  ON  A  0022A 

AUPLIM  =  DESIRED  UPPER  LIMIT  ON  A  0022E 

COMMON  FLAM.TIMF  0023 5 

TYPE  DOUBLE  FLAM  0024  ; 

COMMON  PTRM, GAMMA, PR»PRFMA, ID  0025 

DIMENSION  BC ( 17 ) , A ( 17 ) »FBC(17)  0026 

DIMENSION  PTRM (17, 17, 2) , GAMMA (17,7,17) »PRFMA (7,7)  0027 

DIMENSION  PR (119, 7, 17) »QI (119) ,QJ(119)  '  0028 

DIMENSION  FJPLOT ( 119 ) , JT( 12)  0029 

C  ENTER  DATA  CARDS  HERE  0030  , 


DATA( (QI ( I),I*l,119)a.2, 161.05) , 102 ( .0) ) 
NA  =  1 6 
ALAM=.01 
ID  x  8 
FLAM*3  .0 
PGAM=P= .4 

IYY  =  13421773  .  . 

TIME  =  .125 
ICYCLE»20 
JCYCLE*4 
AL0LIM»4, 


AUPL IM*6. 

C  END  OE  DATA  CARDS 

AL=AL0LIM+1*  $  AU«AUPLIM  +2. 
UNITTxl. 

N=NA+1 
I AMAX=7 
IALASTsO 
D«ELOATF( ID) 

NLA&Tx^EWNxN 

KT»1 

809  IA«KRAN(AL,AU,IYY) 

IF(KT.EQ.l)  113,115 


0031 

0032 

0033 

0034 

0035 

0036 

0037 

0036 

0039 

0040 

0041 


115  T1=-L0GF{ .000000001  +  RANF (-1 ) ) *2. 30258/ALAM 
IF(TJ-TFLC)130.131»132 

130  T2  =  -'  JGF(  .000000001  +  RANF (-1 ) ) *2. 30258/ALAM 
I F (  T1+T2-TFLC )  230.231.1M 

230  T3=-LOGF( .000000001  +  RANF (-1 )) *2. 30258/ALAM 
IF(T1+T2+T3-TFLCJ  331*331.231 

331  NEWN-NLAST-3  5  GO  T0113 

231  NEWN-NLAST-2  $  GO  T0113 
132  NEWN-NLAST  $  GO  T0113 

131  NEWN-NLAST-1 

113  PRINT  8882 » IA.NEWN 

IF(NEWN-IA)  15.13.13  . 

15  IA«NEWN 

13  IF ( IALAST )  11.12.11 
12  CONTINUE 

C  FROM  THIS  NEXT  STATEMENT  TO  NO.  483  IS  CONCERNED  ONLY  WITH  THE  GRAPH 
DO  482  1-1.12 

482  JT ( I )  *8H 

JT ( 1 )*8HE(A/C)  * 

JT ( 3 )=8HSP0TS  = 

JT ( 5 )  *8H  T  * 

JT (7)*8HJ  VS  QJ 
JT ( 8  )*8HVECT0R 
JT(9)=8H  N  * 

JT  ( 1 1 )  ®8H  'A  * 

.  DO  483  1*1.119  <, 

FI  =  I 

483  FJPLOT ( I ) *FI 
IALAST  = IA 

DO  1235  1*1.17 
DO  1235  J=1 » IAMAX 
DO  1235  K-1.17 
1235  GAMMA { I »J,K)*0*0 

C  AT  THIS  PT  THE  LANDING  TRANSITION  PROBABILITIES  ARE  COMPUTED- 
DO  300  IAA*1.  IAMAX 
DO  301  MI =1 . IAMAX 
IFUAA-MI  J31.32.33 

31  PRFMAt IAA.MI )*0. 

.  GO  TO  301 

32  MM1*MI-1 

PRFMAt IAA.MI )«P**MM1 
GO  TO  301 

33  IAM1*IAA-1 
'MM1 *MI-1 

BC(1)*1.0 

PROD*FLOATF( IAA-MI ) 

DO  50  I P*2 .MI 
AIP*FL0ATF( IP-1) 

PROD  *  PROD  +1.0 
50  BCUP)«PROD*BC(  IP-D/AIP 
I  GO* IAA-MI 

PRFMAt IAA.MI )*  ( BC (M I ) * 1 1 .0-P )**( I  GO) )*P**MM1 

381  eaNffNue 

300  CONTINUE 

C  AT  THIS  PT  THE  MAINTENANCE  TRANSITION  PROBABILITIES  ARE  COMPUTED 
DO  100  IT-1.2 


0042 
0043 
0044 
0045 
0046- 
0047 
0048  . 
0049; 
0050: 
0051  : 
0052  ' 
0053  • 
0054  ! 
0055 
0056  : 
0057  i 
0058  ! 
0059  ! 
0060  ' 
0061 
0062  ; 
0063  j 
0064 
0065  : 
0066  * 
0067  i 
0068  1 
0069  • 
0070  j 
0071 
0072 
0073  i 
0074, 
0075  i 
0076' 
0077  ! 
0078 
0079  ; 
0080  1 
0081  ! 
0082  ; 
0083  j 
.  0084  : 
0085  ! 
0086  < 
0087  ! 

0086  j. 

0089 
0090  i 
0091  ' 
0092  ‘ 
0093 
§09*  : 

0095  i 

0096 

0097 


I F ( I T-l ) 25 *25  *26  0098 

25  TAU  =  TIME  0099 

GO  TO  28  0100 

26  TAU  *  UNITT-TIME  0101 

.28  00101  1*1  ,N  0102" . 

DO  102  I J*1 »N  0103 

IF  (I-IJ)  14*199*17  0104 

199  IF(I-ID)  19,19,1999  0105- 

1999  PTRM(I,IJ,IT)»EXPF{-FLAM*TAU*D)  0106 

GO  TO  102  0107 

14  PTRM( I  *  I J* IT ) *0*  0108 

GO  TO  102  .  0109 

19  FJM1*FL0ATFUJ-1)  0110 

PTRM{ I  ,IJ,IT)=EXPF(-FLAM*TAU*FJM1)  0111 

GO  TO  102  0112 

17  IF(I-ID-l)  1*1*2  0113 

1  BC(1)=1.0  0114 

PROD* FLOAT F (I-IJ)  0115 

DO  10  IP  =2 , I J  0116 

AIP=FLOATF( IP-1)  0117 

PROD  *  PROD  +  1.0  0118 

10  BC(IP)  *PROD*BC ( IP-1 ) /A I P  0119 

ELT*EXPF ( -FLAM* TAU )  0120 

PTRM  (  I  *IJ*IT)=BC( IJ)*<1.-ELT)**( I-IJ)*ELT**C IJ-1)  0121 

GO  TO  102  0122 

2  IF ( I J-l-ID)  22,24*24  0123- 

22  CONTINUE  0124 

CALL  PID(I,IJ»IT)  0125 

GO  TO  102  0126  . 

24  D*FLOATF ( ID )  ■  0127 

ELDT«EXPF(-D*FLAM*TAU)  0128 

FACT  *  1.0  0129 

A  C 1 ) *1 .0  0130 

MM  *  I-IJ  0131 

DO  20  M=2 »MM  0132 

FACT*FACT+1 .0  0133 

20  A(M)*A(M-1)*FACT  0134 

201  PTRM( I , IJ»IT)*(D*FLAM*TAU)**< I-I J)*ELDT/A( I-IJ)  0135 

102  CONTINUE  *  0136  . 

101  CONTINUE  0137 

100  CONTINUE  0138 

GO  TO  120  0139 

11  CONTINUE  0140 

IF{  IA-IALA.ST)  120*121*120  0141 

121  IF(NEWN-NLAST) 111*117*111  0142 

117  IALPH* IALASTMl  $  GO  TO  801  0143 

120  CONTINUE  0144  * 

C  AT  THIS  POINT  THE  LAUNCHING  TRANSITION  PROBABILITIES  ARE  COMPUTED  0145 

DO  204  IFF*1»N  0146 

IGM  *  XMINOF  ( IA* IFF)  0147  * 

DO  203  IG*1 ♦ IGM  0148 

IGM1*IG-1  0149 

00  tGB  IH* 1 »N  0156 

I HM1* I H-l  01S1 

BPROD* ( ( 1 .-PGAM ) ** I GM1 ) *( PGAM#* IHM1 )  0152 

86  IF(IG-IA)  91*87,84  0153 


91  IF ( IG+IHM1-IFF)  84,82*84  0154 

87  I F ( IG+IHM1-IFF)85*85*84  0155 

84  GAMMA { IFF*IG*IH)=0.  0156 

GO  TO  202  0157 

82  BC(1)=1.0  0158 

PR0D*FL0ATF( IFF-IG)  0159 

DO  30  I P*2 » IG  0160 

AIP=FLOATF(  IP-1)'  0161 

PROD  »  PRO0  +  1.0  0162 

.30  BC ( I P ) =PROD  *  BC( I  P-1 ) / A I P  0163 

IHM1*IH-1  0164 

TFMP=  PGAM**IHM1  .  0165 

TFMP1=< l.-PGAM)**IGMl  0166 

BPROD  =  TEMP*TEMP1  0167 

GAMMA { IFF  *IG* IH)»BC( IG)*BPROD  0168 

GO  TO  202  -  0169 

85  FBC ( 1 ) *1.0  0170 

PROD=FLOATF ( IGM1-1 )  0171 

DO  40  IP»2* IH  0172 

AIP  =  FLOATFCIP-1)  0173 

PROD  =  PROD  +1.0  0174 

40  FBC ( IP ) *PROD*FBC( I  P-1 1 /AIP  0175 

GAMMAt IFF*IG* IH)*FBC( IH)*BPROD  0176 

202  CONTINUE  0177 

203  CONTINUE  0178 

204  CONTINUE  -  0179 

C  REMOVE  CARDS  FROM  HERE  TO  NO  999  IF  PRINT  OUT  NOT  DESIRED  0180 

PRINT  9»(((I*IJ»IT*PTRM(I*IJ*IT>»IT=1>2)»IJ*1»N)»I»1*N)  0181 

9  FORMAT  ( 1H1/ 1 2 ( 6H  PTRM ( 1 2 » 1H» 1 2 » 1H ♦ 1 2 »3H)  *  E14.5)))  0182 

PRINT  99* ( ( ( IFF*IG»IH*GAMMA( IFF.IG*IH) *IFF*1.N) *IG«1*IA) *IH«1*N>  0183 

99  FORMAT ( 1H1/ ( 2 (7H  GAMMA ( 1 2 » 1H , 1 2 » 1H *  1 2*3H )  •  E14.5)))  0184 

PRINT  999* ( ( I AA*MI »PRFMA( IAA»MI)*IAA*1»I AMAX ) *MI»1 *  I AMAX)  0185 

999  FORMAT  1 1H1/C2C  7H  PRFMA ( 12* 1H* 1 2 *3H )  *  E14.5)))  0186 

C  NOW  THE  TRANSITION  MATRIX  MUST  BF  ZEROED  0187 

111  CONTINUE  .  0188 

DO  899  J31 *  119  0189 

DO  899  K=1 *7  0190 

DO  899  L*l» 17  0191 

899  PR ( J*K* L) *0.0  .  0192 

C  START  COMPUTING  THE  ELEMENTS  OF  EACH  ROW*  I«IX+  (ALPHA  -  1)  X  TOTAL  A/C  0193 
DO  1000  IALPH«1»  IALAST  0194 

801  CONTINUE  0195 

DO  1100  IX*1 »NLAST  0196 

C  COMPUTE  THE  P  ELEMENTS  OF  THE  IAPH*IX  ROW  AND  SUM  THE  ROW  0197 

TSUM*0.  0198 

I=IX+( IALPH-1 )+N  0199 

DO  800  IBETA*1 * IA  ?  0200 

RSUM*0.0  0201 

DO  900  IY*1 *NEWN  0202 

PR( I » I  BETA* IY)*0*  0203 

ILIM*NEWN-I ALPH-IX+2  0204 

PSUM«0.0  0205 

$UM*0.0  0201 

SUML«0.0  0207 

DO  500  IL«1*ILIM  0208 

KLIM«IX+IL-1  0209 


IXPIL  =  IX  +IL  -  1  0210 

SUMM=0.  0211 

DO  600  MI=1,IALPH  0212 

SUMK=0.  0213' 

’  DO  700  IK=1 »KLIM  0214 

IKPMI  =  IK  +MI-1  0215 

IF(IXPIL-NEWN)  701*701*700  0216 

701  I F ( I KPMI-NEWN )  702*702*700  0217 

702  GAMH “GAMMA ( ILIM*IBETA»IL)  0218 

PTRMH1  =  PTRM( IXPIL* IK.l )  0219 

PRFMAH  =  PRFMAt IALPH.MI )  0220 

PTRMH2  “  PTRM ( I KPM I  *  I Y  *2 )  0221 

SUM  =  GAMH  *  PTRMH1  *  PRFMAH  *  PTRMH2  0222 

SUMK“SUMK+SUM  0223 

PSUM=PSUM+SUM  0224 

700  CONTINUE  0225 

SUMM  =  SUMM  +  SUMK  0226 

600  CONTINUE  0227 

SUML  =  SUML  +SUMM  0228 

PSLV.2  “SUML  '  0229 

500  CONTINUE  0230 

RSUM=RSUM+PSUM  0231 

PR( I»IBETA»IY)“PSUM  0232 

900  CONTINUE  0233 

TSUM=T  SUM+RSUM  0234 

800  CONTINUE  0235 

PRINT  888  »  TSUM* IALPH* IX  0236 

888  FORMAT  (  7H  TSUM  *  E15. 8*215)  0237 

1100  CONTINUE  '  0238 

1000  CONTINUE  0239 

C  REMOVE  CARD  FROM  HERE  TO  889  IF  P  MATRIX  PRINT  OUT  NOT  DESIRED  0240 

DO  889  J= 1  *  17  0241 

DO  889  K= 1 ♦ 7  0242 

D0889  L*1 *  17  0243 

I =J+ (K-l ) *N  0244 

889  PRINT  890,(PR(  I,LP*L) ,LP“1*IAMAX) *K,J*L  0245 

890  FORMAT (7E14.5*2HJ=I2»5HK  =  1 *A *2HL“ 1 2 )  0246 

DO  898  1=1*119  0247 

898  QJ( I )“0.0  0248 

C  NOW  MULTIPLY  QI  AND  P  TO  GET  QJ  0249 

805  PRINT  807 *KT* i ALAST ♦ IA  0250 

807  FORMAT ( 1H1.13HQ  VECTOR  CASE  13///  15,15)  0251 

DO  802  IBETA=1 *7  0252 

DO  902  I Y= 1 , 1 7  0253 

CAT  THIS  POINT  CALCULATE  THE  ( IBETA,  IY)TH  ELEMENT  OF  THE  QJ  VECTOR  0254 

J=IY+( IBETA-l )*N  0255 

0P1»0.  0256 

OP»0.  ,  0257 

DO  2001  I ALPH“1 *7  0258 

DO  2201  I X* 1  *  17  0259 

I = I X+( I ALPH-1 ) *N  0260 

QP1 *QI ( I ) #PR ( I  * IBETA , I Y )  0261 

qp=qp»9Pl  0262 

2201  CONTINUE  0263 

2001  CONTINUE  0264 

QJ ( J ) «QP  0265 
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PRINT  8882.IBETA.IY, J.QP  0266 

8882  FORMAT { 2  1 4. 4H  QJCI3.3H  )  =  E14.8)  0267 

90?  CONTINUE  0268 

80?  CONTINUE  0269 

C  CHECK  THE  SUM  OF  THE  0  VECTOR  •  0270 

QSUM=0 •  0271 

00  808  J=  1 . 1 19  0272 

808  QSUM=Q J ( J ) +QSUM  0273 

PRINT  8883 »OSUM  0274 

3883  FORMAT (6H  QSUM=  E15.9)  0275 

DO  333  I  •*  18 » 119  0276 

K  =  (I-D/17  .  0277 

FK=FLOATF(K)  0278 

FMEAN=  FK*QJ ( I ) +FMEAN  '  0279 

333  CONTINUE  0280 

TFLC=FMEAN*FLOATF( JCYCLE)  0281 

PRINT  335 »FMEAN  0282 

335  FORMAT {  17HMEAN  A/C  FLYING  =  F10.4)  0283 

C  STATEMENTS  FROM  THIS  POINT  TO  THE  CALL  DRAW  STATEMENT  REFER  TO  GRAPH  0284 
JT ( 2  )  =  I  CODE ( FMEAN )  0285 

JT  ( 4 )  =  I  CODE ( D )  0286 

FKT=FLOATF(KT)  0287 

JT ( 6 ) =  I  CODE ( FKT )  0288 

FN=FLOATF(NEWN-l )  0289 

JT ( 10 ) = ICODE( FN )  0290 

FIAA=FLOATF( IA-1)  0291 

JT ( 12)=IC0DE(FIAA)  0292 

CALL  DRAW ( 1 19 » F JPLOT  »QJ .0 .0  »4H  ♦ JT *0 *0  * 0 * 0 ,0 .0 ,8  *  8 *0 .LAST  )  0293 

FMEAN  =  0.  0294= 

CNEXT  WE  MUST  MULTIPLY  OJ  AND  P  TO  GpT  QK  AND  SO  ON. • . ( QK+. . .N )  0295 

KT  =  KT+1  0296 

IF (KT-ICYCLE)  803.803*806  0297 

803  DO  804  1=1.119  0298 

804  Q I ( I)=QJ( I )  0299  • 

I ALASTM1=IALAST  0300  i 

IALAST  = IA  0301 

NLAST=NEWN  0302 . 

GO  TO  809  0303  : 

806  STOP  06  0304  ! 

f ND  0305 

SURROUTINF  PID(I.J.IT)  0306' 

COMMON  FLAM, TIME  0307! 

COMMON  PTRM.GAMMA ,PR .PRFMA, I D  0308s 

TYPE  D0U8LF  RC.RDC.PROD  .DID3,DID4,DID5,DEXP  0309  : 

TYPE  DOUBLE  DAN.DID1 .DID2.SUM.DN.ANM1 .FAC. COF.PSUM.PTR. FLAM. TAU.D  0310 
DIMENSION  PTRM( 17.17.2 ) »BC( 11) »BDC( 11 )  0311 

DIMENSION  GAMMAt  17.7.17).  PRFMA(7.7),  PRU19.7.17)  0312 

D=FLOATF ( ID )  0313 

IDP1=ID+1  0314 

IF ( IT-1 )25. 25*26  0315 

25  TAU  =  TIME  0316 

GO  TO  28  0317  ! 

26  TAU=  l.-TIME  0318 

28  CONTINUF  0319 

IMDP1® I-ID  0320 

PTR*0.0  0321 
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PSUM=0. 

DO  200  NJ  s  J»  ID 

DEVELOP  N  TAKEN  J  AT  A  TIME  AND  D  TAKEN  N  AT  A  TIME 

BC(1)=1.0 

PROD=FLOATF(NJ-J) 

DO  10  I P«2  *  J 
ATP*FLOATF( IP-1) 

PROD  =  PROD+1.0 
•10  BC(  IP)=PR0D*  BC(IP-1)/AIP 
BDC { 1 ) *1 *0 

PROD=FLOATF( IDP1-NJ)  . 

DO  20  IQ'=2»NJ 
AIQ=FLOATF( IQ-1) 

PROD=PROD+1.0 

20  BDC( IQ) =PROD*8DC( IQ-1) /AIQ 

COF=BC ( J ) *BDC ( N J ) * ( ~ 1 ) ** ( NJ- J ) 

ANM1=FL0ATF ( NJ-1 ) 

DAN=D/ (D-ANM1 ) 

DTD4=DFXP (-FLAM*TAU*ANM1 ) 

DID1={ DAN**{ I-IDP1 ) )*DID4 
SUM*0. 

DN=0. 

DO  201  K*1 » IMDP1 
FAC=1* 

KM1=K-1 
PROD=0 • 

DO  11  I K=1 ♦ KM1 
PROD=PROD+l • 

11  FAC=FAC*PROD 
IMI DK= I-ID-K 

SUM=( ( FLAM*D*TAU)**KM1 )*DAN**IMIDK  /  FAC 

201  DN=DN+SUM 

DID3=DEXP(-FLAM*D*TAU) 

DID2=DN*DID3 
DID5=DID1-DID2 
PSUM=COF*DID5 
200  PTR  =PTR  +PSUM 
103  CONTINUE 

PTRM( I » J » I T ) *PTR 
102  CONTINUE 
101  CONTINUE 
END 

FUNCTION  KRAN( A*Bt IY ) 

THIS  ROUTINE  RETURNS  AN  UNIFORMLY  DISTRIBUTED  RANDOM  INTEGER 

THIS  ROUTINE  RETURNS  A  INTEGER  RANDOM  NUMBER  .GE,  TO  A 

AL**BOTTOM  LIMIT  (INCLUDED)  FOR  THE  RANDOM  NUMBER 

p  ,  T0P  LIMIT  (NOT  INCLUDED)  FOR  THE  RANDOM  NUMBER 

SET  IY  ONLY  ONCE  IN  MAIN  PROGRAM  FOR  EACH  SET  OF  RANDOM  NUMBERS 

9QMI  GOOD  STARTING  VALUES  FOR  IY  FOLLOW 

13421773 

33554433 

8426219 

42758321 


03  22  \ 

0323 

0324 

0325  : 

0326  * 

0327  , 

0328 

0329  * 

0330 

0331 

0332 

0333 

0334  • 

0335 

0336 

0337  • 

0338  • 

0339 

0340 

0341  , 

0342  . 

0343  . 

0344  ; 

0345  j 

0346  i 

0347  .. 

0348  ; 

0349 

0350- 

0351 

0352 

0353 

0354 

0355 

0356 

0357 

0358 

0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0368 

0369 

0370 

0371 

0372 

0373 

0974 

0375 

0376 

0377; 
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uuuuuu 


rtrwv  -.-Mi-la 


56237485  0378 

62104023  0379 

ANY  OF  THESE  MAY  BE  USED  •  0380 

0381 

THIS  ROUTINE  MAY  BE  USED  IN  FORTRAN  60  OR  63  0382 

,  0383 

IY  «  3125  *  IY  0384 

IY  *  IY  -  UY/67108864)  *  67108864  0385 

FY  »  IY  0386 

-  KRAN  =  FY/67108864.  *  (B-A)  +  A  0387 

RFTURN  '  0388 

END  .  0389 

FINIS  .  0390 

-EXECUTER.  0391 
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0  VECTOR  CASE 


1 


5  5 

1  1  OJ(  l  )  =  2. 647 62 77:>E-o2 
1  2  OJ(  2  )=3. 020208316-02 
1  3  QJC  3  >=1.724368230-02 
1  4  QJ(  4  >=6.681525100-03 
1  5  QJ(  5  >=2.02525670E-03 
1  6  0 J (  6  >=5.297568230-04 
1  7  QJ(  7  >=1.300020330-04 
1  8  0 J (  8  )=' .22530060E-05 
1  9  QJ<  9  )  =  8 . 588570 0  QE-06 
1  10  QJC  10  >=2.862856676-06 
1  11  OJ<  11  >=8.349998610-07 
1  12  QJC  12  >=2. 087499650-07 
1  13  QJC  13  >=4.348957610-08 
1  14  QJ(  14  >=7.248262686-09 
1  15  QJ(  15  >=9.060328360*10 
1  16  QJC  16  >=7.550273630-11 

1  17  Q J (  17  >=3.145947350-12 

2  1  OJC  18  >=3.025194336-02 
2  2  QJC  19  > =3. o738l 5526-02 
2  3  OJ(  20  >=1.543854056-02 
2  4  Qj(  21  >=5.188651076-03 
2  5  Q j (  22  >=1.343744376-03 
2  6  QJ(  23  >=2. 961643020-04 
2  7  QJ(  24  >=6.072270210-05 
2  8  QJC  25  >=1.262036490-05 
2  9  QJ<  26  >=2.862895036-06 
2  10  0 J <  27  >=8.35oiio5l0-o7 
2  11  0 J (  28  >=2.087527636-07 
C  12  OJC  29  >  =4 . 349015896-08 
•  33  QJ(  30  >=7.248359816-09 

**  OJ(  31  >=9.060449776-10 
2  15  QjC  32  >=7.550374810-11 
2  16  0 J (  33  >=3.1459895QE-12 

2  "  17  Q J (  34  >=  0 

3  1  0 J (  35  >=3,409690436-02 
3  2  QJ(  36  >=3.075975310-02 
3  3  QJ(  37  >=1.352651076-02 
3  4  QJ<  38  >=3. 915317810-03 
3  5  QJ<  39  )=8.57o928oiE-04 
3  6  QJ(  4o  )  =  1. 56625 190E*04 
3  7  OJ(  4i  )=2.6227io346-o5 
3  8  OJ(  42  >=4.44ip5l686-o6 
3  9  QJ(  43  >=8.351259556-07 
3  10  QJC  44  >=2.0878l489E-07 
3  11  GJ(  45  >=4.349614386-08 
3  12  QJ(  46  >=7.249357256-09 
3  13  0 J (  47  >=9.061696560-10 
3  1<*  QJC  ia  >  =  7.55l4i380E-ll 
3  15  QJC  49  >=3.145422426-12 
3  16  QJ(  50  >=  0 
3  17  0 J (  51  )=  0 


4 

1 

QJ( 

52 

>=3.797642956-02 

4 

2 

QJ< 

53 

>=3.032240616-02 

4 

3 

QJC 

54 

>=1.1623o913E-02 

4 

4 

QJC 

55 

)  =2 . 87808733E-Q3 

4 

5 

QJC 

56 

)=5.26888829E-04 

4 

6 

QJC 

57 

>=7.847204546-05 

4 

7 

QJC 

58 

>=1.045l6559E-05 

4 

8 

QJC 

59 

)=1.39l45026E-06 

4 

9 

QJC 

6o 

>=2.089606240-07 

4 

10 

QJC 

61 

>=4.353346346-08 

4 

11 

QJC 

62 

>=7.255577230-09 

4 

12 

QJC 

63 

>=9.069471540-10 

4 

13 

QJC 

64 

>  =  7.557892956-h 

4 

14 

QJC 

65 

>=3.149122060-12 

4 

15 

QJC 

66 

>=  o 

4 

16 

QJC 

67 

>=  o 

4 

17 

QJC 

68 

>=  o 

5 

1 

QJC 

69 

>=4.916428330-01 

5 

2 

QJC 

70 

>=1. 437181270-01 

5 

3 

QJC 

71 

>=2.695318530-02 

5 

4 

QJC 

72 

>=3.787099626-03 

5 

5 

QJC 

73 

>=4.237326536-04 

5 

6 

QJC 

74 

>=3.9938l550E~05 

5 

7 

QJC 

75 

>=3.419140620-06 

5 

8 

QJC 

76 

>=2.985611350-07 

5 

9 

QJC 

77 

>  =  3 . 11088104E-08 

5 

10 

QJC 

78 

>=4.951147550-09 

5 

11 

QJC 

79 

)  =5 . 9210  8046E-10 

5 

12 

QJC 

80 

>=4.728868746-11 

5 

13 

QJC 

81 

)=1.8913981QE-12 

5 

14 

QJC 

82 

>=  a 

5 

15 

OJC 

83 

>=  o 

5 

16 

QJC 

84 

)=  o 

5 

17 

QJC 

85 

')=  o 

6 

1 

QJC 

86 

)=  o 

6 

2 

QJC 

87 

>=  o 

6 

3 

QJC 

88 

>=  0 

6 

4 

QJC 

89 

>=  o 

6 

5 

QJC 

90 

>=  o 

6 

6 

QJC 

91 

)=  o 

6 

7 

QJC 

92 

>=  o' 

6 

8 

QJC 

93 

)=  o 

6 

9 

QJC 

94 

>=  o 

6 

10 

OJC 

95 

>=  o 

6 

11 

QJC 

96 

>=  o 

6 

12 

QJC 

97 

>=  o 

6 

13 

QJC 

98 

)=  o 

6 

14 

QJC 

99 

)=  o 

6 

15 

OJC100 

>=  o 

6 

16 

QJC 101 

>=  0 

6 

17 

QJC102 

)=  o 
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7 

1 

QJ(103 

)  = 

0 

7 

2 

QJCL04 

)  = 

0 

7 

3 

QJ(105 

)  = 

0 

7 

4 

Q J ( 106 

>  = 

0 

7 

5 

Q J  < 1 07 

)s 

0 

7 

6 

QJ( 108 

>  -  • 

0 

7 

7 

0 J \ 109 

)S 

0 

7 

8 

0.(110 

)  = 

0 

7 

9 

QJdii 

)  = 

0 

7 

10 

0 J ( 112 

)  = 

0 

7 

11 

Q J ( 113 

)  = 

0 

7 

12 

0J(H4 

)=  . 

0 

7 

13 

OJ(ll5 

)  = 

0 

7 

14 

0J(ll6 

)  = 

0 

7 

15 

0 J ( 117 

>  = 

0 

7 

16 

0 J ( 118 

)  = 

0 

7 

17 

0 J ( 119 

)  = 

0 

QSUM 

*1.1 

DOOOOOOOOE  00 

EAN 

A/C 

flying 

s3.i666E  00 

gkaph  TITLED 

E ( A/C )  s3.17E^OOSPOTS  =  8.00E--0Q 
J  VS  OJ  VECIOK  N  s  1 . 608*01 

HAS  BEEN  PLOTTED. 


51 


T  =  1 . OOE+OO 
A  s  4. 0OE*OO 


v  i 


APPENDIX  IV 


SAMPLE  RESULTS 

The  following  pages  present  the  values  of  the  elements  of  *he  proba¬ 
bility  distribution  vector  (Q  J)  and  its  graphical  plot  for  five  consecutive 
iterations,  i.  e. ,  Q  x  Pn  for  n  =  1,  2,  . . .  ,  5.  The  inputs  are  those 
shown  on  the  first  page  of  Appendix  III  between  statement  No.  30  and 
No.  31.  The  printouts  of  the  transition  matrices  and  their  computational 
elements  are  omitted.  The  plot  was  made  using  the  DRAW  subroutine  in 
the  U.  S.  Naval  Postgraduate  School  computer  facility  library.  Each 
vector  printout  contains  the  values  of  all  119  states  possible  (7  x  17)  and 
is  headed  by  the  past  value  of  A  +  1  and  the  next  value  of  A  +  1 .  The  two 
indices  preceding  each  element  represent  fj  +  1  and  j  +  1,  in  the  notation 
of  section  3.  For  example,  in  the  first  row  on  the  next  page,  the  "1  1" 
indicates  that  the  probability  of  being  in  state  (0,  0)  after  one  iteration 
is  =  .  026,  where  the  value  of  A  is  4  over  the  first  iteration.  Each  graph 
is  labeled  with  the  expected  value  of  a/c  flying,  the  number  of  mainte¬ 
nance  spots  available,  the  vector  number  (T),  total  number  of  a/c  avail¬ 
able  (N),  and  the  desired  number  of  a/c  on  station  (A).  The  "E"  notation 
indicates  the  power  of  10  to  multiply  by.  This  sample  run  demonstrates 
the  loss  in  total  a/c  and  variable  a/c  on  station. 
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002 


H*SCALE  -  2.00E+81  UNITS-' INCH. 

Y -SCALE  -  1J30E-W  UNITS'INCH. 

E(fVC)  =3. .17E+00SPOT3  =  8.00E+00  T  = 
J  US  QJ  UECTOR  N  -  1.60E+01  A 


1.00E 

■  4lo0 


0  VECTOR  CASE  2 


5  5 

1  1  QJ(  1  )=9.12184624E-Q6 

1  2  QJ(  2  >=7.944910286-06 

1  3  QJ(  3  >=3.413805336-06 

1  4  Q J (  4  )=9.839l3729E-07 

1  5  0 J (  5  >=2.205303596-07 

1  -6  QJ(  6  ) =4 . 28567319E- 08 

1  7  0 J (  7  )=7.97994342E-09 

1  8  0 J (  8  )=1.56348438E-09 

1  9  0 J (  9  )=3.46529220E-10 

1  10  QJ<  10  >=1.012212426-10 

1  11  0 J (  11  >=2.581730036-11 

1  12  0 J <  12  )=5-.633852016-l2 

1  13  OJ(  13  )=1. 022956166-12 

1  14  QJ(  14  )=1.48426307E-13 

1  15  QJ(  15  )=1.61419113E-14 

1  16  QJ<  16  )=1.17o3o587E-15 

1  17  QJ(  17  )s4.24536876e-17 

2  1  0 J (  18  )=1.69lo5396E-o4 

2  2  QJ(  19  )=1.3io58642E-0^ 

2  3  0 J (  2o  >=4.93425o7iE-o5 

2  4  QJ(  21  )=1.224i68l9E-05 

2  5  OJ(  22  >=2.3i695297e-()6 

2  6  0J(  23  )=3.73462l97E-o7 

2  7  0 J (  24  )=5.7o66i929E-08 

2  8  0J(  25  )=9.22l54646E-09 

2  9  OJ(  26  J-1.72702939E-09 

2  10  OJ(  27  )=4.4l48l68lE-10 

2  11  QJ(  28  >=9.658696766-11 

2  12  OJ(  29  >=1.75883018E-11 

2  13  OJ(  30  )=2.56o23495E-12 

2  14  Q J (  31  )  =2 . 794340016-13 

2  15  0 J (  32  )=2. 03393035E-14 

2  16  Q J (  33  )=7.41007441E-16 

2  17  OJ(  34  >=  0 

3  1  QJ<  35  >=1.425253196-03 

3  2  0 J (  36  )=9.83757697e-04 

3  3  0 J(  37  )=3.2433l35uE-04 

3  4  OJ(  38  >=6. 902340106-05 

3  5  0J(  39  )=1.09403657E-05 

3  6  QJ(  4o  )  =1 . 44o32446£-06 

3  7  QJ<  4i  )  =  1 .76353392E-07 

3  8  QJ<  42  )=2.28l3l763E-08 

3  9  0 J (  43  )=3. 510782906*09 

3  10  QJ(  44  )s7.7o559333e-10 

3  11  QJ(  45  )=1.4085ll7l£-lo 

3  12  0 J (  46  >  =  2.'o5928494E-ll 

J  3  13  Q J (  47  )  =  2. 258692916-12 

'  3  14  Q J (  48  >=1.65302266E-13 

3  15  QJ<  49  )s6. 057991 94E-15 

3  16  0 J (  50  >=  0 

3  17  Oj(  5l  )=  0 


4  1  OJC  52  )  =  7. 230352426-03 
4  2  OJC  53  )=4.45693238E-03 
4  3  0 J (  54  )si. 289666766-03 
4  4  0 J (  55  )=2. 355005326-04' 
4  5  OJ(  56  )=3. 111562326-05 
4  6  Q J (  57  )=3.30021596E-06 
4  7  QJ<  58  )=3.l5l23l46E-07 
4  8  QJ(  59  )s3.l3674353E-08 
4  9  QJ(  6o  )=3.810^4l7lE-09 
4  lo  QJ<  6l  )=6. 986635916-10 
4  11  QJ<  62  )=1. 025687386^10 
4  12  QJ<  63  )  =  1.i3o922556-H 
4  13  QJ<  64  )=8.32962386E-1 3 
4  14  QJ<  65  )=3.o755o593E-l4 
4  15  OJ(  66  )=  0 
4  16  Q J (  67  )s  0 

4  17  QJ(  68  )=  0 

5  1  Q J (  69  )s7.70076632E-01 
5  2  QJ<  7o  )=1.88l59529E-01 
5  3  QJ<  7l  )=2. 324804126-02 
5  4  QJ<  72  )=1.93788734E-o3 
5  5  QJ<  73  )=i. 227620686-04 
5  6  Q J (  74  )=6. 394207946-06 
5  7  Q J (  75  )=3.02738808E-07 
5  8  QJ<  76  >=1.563l2395>E-08 
5  9  QJ<  77  >=1.135192056-09 
5  10  QJ<  78  >=1.50627061E-10 
5  11  QJ<  79  )=1.51200^91E-11 
5  12  QJ<  80  >=1.020934346-12 
5  13  OJC  81  >=3.478808656-14 


5  14  Q J (  82  >=  0 
5  15  QJ<  83  >=  0 
5  16  QJ<  84  >=  0 

5  17  OJC  85  )=  0 

6  1  QJ<  86  )=  0 
6  2  QJ(  87  )=  0 
6  3  QJ<  88  )=  0 
6  Qj(  89  )=  0 
6  5  OJC  9o  >=  0 
6  6  QJ<  91  )=  0 
6  7  QJ<  92  >=  0 
6  8  OJC  93  );  0 
6  9  QJ<  94  )=  0 
6  10  QJ(  95  )s  0 
6  11  OJC  96  )=  0 
6  12  QJ<  97  )=  0 
6  13  QJ(  98  >=  0 
6  14  QJ(  99  >=  0 
6  15  QJ<100  >=  0 
6  16  QJ(101  >=  0 
6  17  QJ( 102  >=  0 
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7  1  G)Jd03  )=  0 

7  2  QJ(104  )=  0 

7  3  QJ( 105  )s  0 

7  4  QJ(106  )=  0 

7  5  QJ( 107  )=  0 

7  6  0 J ( 108  )=  0 

7  7  Q J ( 109  )=  0 

7  8  OJ(110  )*  0 

7  9  OJ(lll  )=  0 

7  10  QJdl2  )=  0 

7  11  0 J ( 113  )=  0 

7  12  QJdl4  )=  '  0 

7  13  QJdl5  )=  0 

7  14  QJdl6  )-  0 

7  15  0 J (117  >=  0 

7  16  QJ(118  )  =  0 

7  17  QJ(119  )=  0 
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000  002  005  006 


006  010 


/ 

X -SCALE  -  2.00Ev01  UNITS'INCH.  { 

Y -SCALE  -  1.00E-01  UNITS'INCH. 


E(A/C)  =3»98E+00SPOTS  =  8.00E+00 
J  US  QJ  UECTOR  N  =  1.60E+01 


T  »  2.00E 
A  =  1.01 
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0  VECTOR  CASE  3 


5  7 

1  1  OJ(  1  >--2.835683l4fci~{p 

1  2  0 J (  2  )=2.l585i38iE-0,j 

1  3  0 J (  3  )  =  7 . 9^343376e-(j6 

1  4  QJ(  4  >=1.9l420746E-o* 

1  5  OJ(  5  )=3.49l25464g-07 

1  6  QJ<  6  )=5.37405i66e-08 

1  7  QJ(  7  )=7.77939446E-09 

1  8  OJ(  8  )  =  1 . l87566o6E-0 9 

1  9  QJ(  9  ):2.1128384lg-i0  , 

1  10  OJ(  10  >=5.i899266uE-ll  ' 

1  11  QJ<  11  )=1.0852862UE-11 

1  12  OJ(  12  )  =  1'.87739309E-12 

1  13  0 J (  13  )=2.57762863E-13 

1  14  QJ<  14  >  =  2 . 631762916-14 

1  15  QJ<  15  )=1.77499822E-15 

1  16  QJ<  16  )  =  5 . 92695364E-17 

1  17  QJ<  17  )  =  0 

2  1  QJ(  18  )=4.8l7ll691E-04 

2  2  QJ(  19  >=3.26i73533E-04 

2  3  OJ(  20  >=1.04905695E-04 

2  4  Q J (  2i  )-2.l6i42330E-0* 

2  5  OJ(  22  )=3.28306882E-06 

2  6  OJ(  23  >s4.o9o33643E-Q7 

2  7  OJ(  24  >=4.68o83562E-03 

2  8  QJ(  25  >=5.625863736-09 

2  9  OJ(  26  >=8.102581696-10 

2  10  OJ(  27  >=1. 6931773 6E -10 

2  11  Q J (  28  )=2.92764116E-11 

2  12  0 J C  29  )  =  4.0l87i703E-l2 

2  13  OJ(  3o  )=4.i0297009E-l3 

2  14  GJ<  3l  )-2.76746676E-l4 

2  15  OJ(  32  T-9.24174352E-16 

2  16  OJ(  33  >  =  0 

2  17  OJ(  34  >  =  0 

3  1  QJ(  35  )  =3 .675899586-03 

3  2  Q J (  36  )=2.220293l3E-03 

3  3  0 J (  37  )=6.25702594E-04 

3  4  QJ(  38  >=1 .10292585E-04 

3  5  QJ(  39  >  =  1 . 3890038 7E-05 

3  6  OJ(  4o  )=1.38024736E-06 

3  7  Q J (  41  ) =1 .210573586-07 

3  8  OJ(  42  >=1.092620206-08 

3  9  0 J <  43  )-l. 2136293 7E-09 

3  10  OJ(  44  )=2. 090576196-10 

3  11  Q J (  45  )=2.86o5573iE-ll 

3  12  QJ(  46  >=2. 913667126-12 

3  13  QJ(  47  )=1. 962826196-13 

3  14  GJ<  48  )s6. 555493746-15 

3  15  0 J (  49  )=  0 

3  16  0 J (  5o  )s  a 

3  17  QJ(  5i  )s  0 


4  l  GJ<  52  )si. 664993696-02 
4  2  QJ(  53  )=9.0l02324VE-03 
4  3  QJ(  54  >=2.238421466-03 
4  4  QJ<  55  >=3.400063736-04 
4  5  QJ(  56  >=3.571409956-05 
4  6  OJ(  57  >=2.825935626-06 
4  7  QJ(  58  >=1.858944676-07 
4  8  QJ<  59  >  =  1 . 19355503E-08 
4  9  QJ(  60  >=9.515295576-10 
4  10  OJ(  61  >=1.285381466-10 
4  11  QJ(  62  >=1 .28940846E-11 
4  12  Q J (  63  ) =8 .533209246-13 
4  13  0 J (  64  )=2.792626?0E-14 


4  14  QJ(  65  >=  0 
4  15  QJ<  66  )=  0 
4  16  OJ(  67  )=  0 
4  17  QJ<  68  >=  0 


5  1  QJ<  69  >=4.982535976-02 
5  2  QJ<  7o  >=2.427388326-02 
5  3  QJ(  71  >=5.358957036-0? 
5  4  QJ<  72  >=7.lc2829736-04 
5  5  QJ<  73  )=6.3358853oe-03 
5  6  QJ<  74  >=4.082328686-06 
5  7  QJ<  75  >=2.o5o24689E-q7 
5  8  QJ<  76  >=9.267590696-09 
5  9  QJ<  77  )=5.o446l339E-10 
5  10  QJ<  78  >=5.l326l492E-ll 
5  11  QJ<  79  >=3.44456222E-12 
5  12  QJ<  80  >=1.1427505dE-13 


5  13  QJ<  81  >=  0 
5  14  QJ<  82  >=  0 
5  15  QJ<  83  >=  0 
5  16  QJ(  84  )=  0 
5  17  QJ(  85  >=  0 


6  1  QJ<  86  >=1.04300023E-01 
6  2  QJ<  87  >=4.528950106-02 
6  3  QJ(  88  >=8.794537716-03 
6  4  OJ(  89  ) =1 . 00690538E-03 
6  5  OJ(  90  >=7.555l2291E-05 
6  6  0 J (  91  >=3.924185766-06 
6  7  G J (  92  >=1.476993196-07 
6  8  QJ<  93  >=4.436304486-09 
6  9  QJ<  94  >=1.438999356-10 
6  10  QJ<  95  >=9.803964676-12 
6  11  GJ<  96  )=3.3oi44o45e-l3 


6  12  QJ<  97  )=  0 
6  13  QJ<  98  )=  0 
6  14  QJ<  99  >=  0 
6  15  OJ(100  >-  0 
6  16  OJ(101  >=  0 
6  17  QJd02  >=  0 
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7 

1 

QJ(103 

)  =  5 . 5563984QE-01 

7 

2 

QJ(104 

)  =  1 .48858751E-01 

7 

3 

QJ<105 

)=1.84227232E-02 

7 

4 

QJ(106 

)=1.37794578E-03 

7 

5 

Q J ( 107 

)  =  6 . 83380  035E “05 

7 

6 

0 J ( 108 

)  =  2 . 32262953E-06 

7 

7 

0 J ( 109 

)  =5 , 43332378E-08 

7 

8 

OJ(110 

)  =  8 . 8065396  7E-10 

7 

9 

QJdll 

)  =  1 . 18536160E-11 

7 

10 

QJdl2 

)=3.82356527E-13 

7 

11 

0J(113 

)=  0 

7 

12 

0 J ( 114 

>=  '  0 

7 

13 

QJ(115 

)=  0 

7 

14 

QJdl6 

)=  0 

.7 

15 

Q J (117 

)=  0 
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)=  0 
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Q  VECTOR  CASE 


4 


7  5' 


1 

1 

Q  J  ( 

1  > 

=4.37066460E-q4 

4 

1 

QJ< 

52 

)=8. 03893389E-02 

1 

2 

QJ< 

2  ) 

=2.66337542E-o4 

4 

2 

QJ< 

53 

)=3.58842827E-02 

1 

3 

QJ( 

3  ) 

=7 . 6i707863E-q5 

4 

3 

QJ< 

54 

)  =  7 . 24458  723E-0  3 

1 

4 

QJ< 

4  ) 

=1.37453763£-05 

4  . 

4 

QJ( 

55 

)=8.757565l3E-04 

1 

5 

QJ( 

5  ) 

=1.796i7995E-06 

4 

5 

QJ( 

56 

) =7 . 09643835E-03 

1 

6 

QJ( 

6  ) 

=1.39ii2432E-07 

4 

6 

0  J  ( 

57 

) =4 . 12757545E-Q6 

1 

7 

QJ( 

7  ) 

=1.8l09889lE-o8 

4 

7 

QJ< 

58 

)=1.85478755E-07 

1 

6 

0  J  < 

8  ) 

=1.843267206-09 

4 

8 

QJ< 

59 

)s7.45697706E-09 
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2uhp reject  ansAer,  system  number*#  task  somber,  eta. 

ft  ORfOtNATOR'8  REPORT  IBHMEE(8):  Enter  the  efft- 
einl  repert  osmber  by  wtdeh  the  deeumemt  wUl  be  Identified 
SOd  sdstrsUed  by  the  crigleetlag  activity.  Tide  amdtr  most 
he  «dpw  te  this  repert* 

#6.  OTHER  REPORT  NUMBER'S):  If  the  repert  hee  heso 
assigned  any  ether  repert  number#  farther  6 f  the  eripmsr 
er  6 y  the  sponsor),  sloo  sotsr  this  uamberis). 

10.  AV  AIL  AKLIT  Y/LHOT  ATTON  NOTICE*  Eotor  aay  lim* 
iutiena  on  further  disosmUstion  of  the  report#  other  then  thee 


impoaod  by  security  classification#  using  standard  statements 
eueh  am 


(t)  “Qualified  regpaeelers  my  ebteio  capias  ef  this 
repert  hem  DDC” 


(3)  “Fereign  eeneuecese 

repert  by  DDC  is  net 


dieeemlantles  ef  tbie 


(3)  “U  8.  Qovsrnmsnt  agon  aloe  an  abtaln  sepias  of 
thU  report  directly  hem  DDC,  Other  gsMMad  PPC 


(4)  <#U.  X  sdiitary  ogsoclss  may  ehtein  espies  of  this 
repert  directly  hem  DDC  Other  testified  gears 
shell  request  through 


(9)  “All  distribution  ef  this  repert  Is  osntieUed.  Qual¬ 
ified  DDC  users  shell  request  through 


cute  this  feet 
1L  WPPLM 


newest  ef  CeemaereA  hr  sale  te  tbs  pdhlie#  imdb 

m  mMmwm 


uoffUTAT  worm  v—tmimtm iwpi 


UWOHwwfi  lnutifiiii, 

13.  AaVTUACf:  fetor  «*  *****  #*■«  «  Mrf  art  fertwl 
tummaryef  fhe  da  cent  at  iodissttss  ef  the  repert#  seen  theatfi 
it  may  ale#  appear  elsewhere  la  ht  bedy  a t  the  teeholcel  re¬ 
port.  If  additional  space  is  regatta  d,  s  esatlsaitien  sheet  shall* 
he  attached. 

It  is  highly  desirable  that  the  ebe tract  ef  classified  reports 
be  uncleeslned.  Each  pore  graph  of  the  ebe  tract  ehell  end  with 
an  indication  of  the  sdiitary  security  cUselAcatiea  of  the  in¬ 
formation  in  the  pampaph#  represented  as  (T9),  (9),  (C),  er  (V). 

There  I*  so  limits  ties  eo  the  long*  ef  the  obstruct  How¬ 
ever,  the  eaflwssted  leaglh  is  hem  190  te  339  wards. 


14.  EEY  IfORDOt  Eey  wards  are  tsshrdsolly  onafcaghit  terms 
or  short  phrases  that  characterise  a  report  sod  map  be  used  sa 
laden  aotriee  far  aafoleglgg  tha  repert  Eey  wards  must  be 
selsatsd  se  hit  ae  security  elaaaiRcatieo  is  resulted.  Idauti* 
Rare,  each  aa  egalpmeot  model  daslapatira#  trade  name,  milHacy 
prafeet  aada  name,  aar^aabio  taaadao#  may  ba  Mad  aa  bey 
wards  bat  will  be  Meeritt  so  ladlsstlsu  ef  tedhaleel  con¬ 
test.  The  asslgamaat  ef  Uahs#  refea#  sod  wei^re  is  aptisorl. 


